SPLITTING METHODS FOR THE NONLOCAL FOWLER EQUATION 



AFAF BOUHARGUANE AND REMI CARLES 



Abstract. We consider a nonlocal scalar conservation law proposed by Andrew C. 
Fowler to describe the dynamics of dunes, and we develop a numerical procedure based on 
splitting methods to approximate its solutions. We begin by proving the convergence of the 
well-known Lie formula, which is an approximation of the exact solution of order one in 
time. We next use the split-step Fourier method to approximate the continuous problem us- 
ing the fast Fourier transform and the finite difference method. Our numerical experiments 
confiiTn the theoretical results. 



1. Introduction 
We consider the Fowler equation ffTTl fT2l : 



^^^^ dtu{t,x) + d^,i^ — jit,x)+I[u{t,-)]{x)^d'Mt,^)^0, xeR,t>0, 

u{0,x) — uq{x), X e R, 

where u = u{t,x) represents the dune height and I is a nonlocal operator defined as 
follows: for any Schwartz function ip E 5(R) and any x e R, 

(1.2) := / \^\-h^"(x-OdC 



Jo 

We refer to |[T]|2l|6l for theoretical results on this equation. 

Remark 1.1. The nonlocal term I is anti -diffusive. Indeed, it has been proved in [jj that 

(1.3) ^™)(e) = -4.^r0) (^i-^sgn(0^^ IC^, 

where denotes the Fourier transform normaUzed in ( |l.ll| i. Thus, I can been seen as a 
fractional power of order 2/3 of the Laplacian, with the "bad" sign. It will be clear from 
the analysis below that our results can easily be extended to the case where I is replaced 
with a Fourier multiplier homogeneous of degree A g]0,2[, as in ||4l, and not only A = 4/3. 

We assume that the initial data uo belongs to i7^(R), and thus ( |l.l| i has a unique solu- 
tion belonging to C([0, t],H^{'R)) for all t > 0, from 1 1 ]. We will denote u{t, •) by 5*un; 
5* maps H^{H) to itself. Duhamel's formula for the continuous problem ( |l.l| l reads 

1 f* 

(1.4) u{t, ■) S*Uo = K{t, .j^uo^-J^ d,K{t - s, •) * {S'u^f ds, 
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where K{t, ■) = ^ (e ''^'^j is the kernel of the operator X — d^, and tpi is defined by 

(1.5) MO - 47r2C2 _ ^^|^|4/3 ^ ,5^^|^|i/3^ 

where ax, hi are positive constants. 

Recently, to solve the Fowler equation some numerical experiments have been per- 
formed using mainly finite difference approximation schemes ||3] ID. However, these 
schemes are not effective because if we opt for an explicit scheme, numerical stability re- 
quires that the time step At is limited by O(Ax^). And, if we choose an implicit scheme, 
we have to solve a large system which is a computationally expensive operation. Thus, 
the splitting method becomes an interesting alternative to solve the Fowler model. To our 
knowledge, there is no convergence result in the Uterature for the splitting method associ- 
ated to the Fowler equation. This method is more commonly used to split different physical 
terms, such as reaction and diffusion terms, see for instance [18]. Splitting methods have 
also been employed for solving a wide range of nonlinear wave equations. The basic idea 
of this method is to decompose the original problem into sub-problems and then to approx- 
imate the solution of the original problem by solving successively the sub-problems. Vari- 
ous versions of this method have been developed for the nonlinear Schrodinger, Korteweg- 
de-Vries and modified Korteweg-de-Vries equations, see for instance lfT3l [TTl [T9l |2T1 . 
For the Fowler model we consider, separately, the linear Cauchy problem 

(1.6) —^I[v{t,-)]-i^dlv^Q- v{{),x) = vo{x), 
and the nonlinear Cauchy problem 

(1.7) — -£C>> = 0; w{0,x)=wo{x), 



where e, rj are fixed positive parameters such that e + rj = 1. Equation ( |1.7| l is simply the 
viscous Burgers' equation. We denote by X* and F*, respectively, the evolution operator 
associated with ( |1.6| i and ( |1.7| i: 

v{t,-) ■.^X*vo = D{t,-)*vo, 

where D{t, •) = (e"*"^^) with = 'in'^vf - ai|f l*/^ _^ bx^\^\^^^, and 

1 /■* 

(1.8) w{t,-):^Y'wo^G{t,-)*wo~-J d^G{t - s, ■) * {Y'wnf ds, 
where G is the heat kernel defined by 

Git, .) = J-i(e-*(4-^-l-l')) = ^^e-^. 

Furthermore, the following i'^-estimate holds 

(1.9) \\Y'w\\l2^j,^ < \H\lhr)- 

Let us explain the choice of this decomposition. First, we can remark that if we do not 
consider the nonlinear term in ( |1.1[ ), the analytical solutions are available using the Fourier 
transform. Thus, the linear part may be computed efficiently using a fast Fourier transform 
(hereafter FFT) algorithms. Note also that the Laplacian and the fractional term I cannot 
be treated separately. Indeed, the equation ut + I[u] = is ill-posed. We next decide to 
handle the nonlinear term by adding a bit of viscosity in order to avoid shock problems in 
the standard Burgers' equation. Therefore, the splitting approach presented in this article 
differs from e.g. the one analyzed in |10|, which corresponds to assuming e = in the 
above definitions. The splitting operators associated to this approach when 1 = (which 
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amount to considering alternatively the heat equation and the Burgers equation) have been 
studied in jTSll (as well as other equations involving the Burgers nonlinearity, such as the 
KdV equation, see also |T4l). The main difference with the result presented in fTSl is 
that the operator X is not a differential operator, so its action on nonlinear terms is rather 
involved , while this point would be needed to compute the Lie commutator between the 
two operators 

A(v) = ^^lv~X\v\, B{v) = -vd^v + e^dlv, 

as in e.g. 031 [TtII . Note that since we consider the Lie splitting operator, we could in 
principle use the exact formula established in |7| for the local error But then again, we 
face the problem to compute [A, B], which is the only term that we cannot estimate in fT] 
Theorem 1]. We finally point out that we use smoothing effects associated to the viscous 



Burgers equation; see Corollaries 3.8 and 3.9 



We motivate this choice by the presence of artificial diffusion in classical numerical 
schemes used to solve the convection equations. An alternative to reduce this effect is to 
consider numerical schemes of high order which are usually computationally expensive 
and do not seem to be very useful for the Fowler model because of the diffusion term. 

We consider the Lie formula defined by 
(1.10) Zi=X^YK 

The alternative definition Z£ ~ Y^X^ could be studied as well, leading to a similar result. 
Also, the following evolution operators 

Zl = or Zl = yV2;^«yt/2^ 

corresponding to the Strang method |20| could be considered. Following the computations 
detailed in the present paper for the case ( |1.10[ ), it would be possible to show that the other 
Lie formula generates a scheme of order one, and to prove that the Strang method is of 
order two (for smooth initial data), in the same fashion as in, e.g., |l5][l7l. This fact is 
simply illustrated numerically in Section|6] to avoid a lengthy presentation. With Z\^ given 
by ( |1.10[ ), our main result is: 

Theorem 1.2. For all uq G (R) and for all T > 0, there exist positive constants Ci , C2 
and Ato such that for all At g]0, Ato] and for all n Cz IS! such that ^ nAt ^ T, 

||(Zf )"uo - 5"^*uo||l2(R) ci At and || (Zf )"uolk3(R) ^ C2. 

Here, ci,C2 and Ato depend only on T, p — max^gjo T] ||'5'*Mo||ff2(j^-|, and \\uq\\ (^ny 

Remark 1.3. It will follow from Lemma p.lll that 

p= max 1 15*^0 11^2 (R) ^ Ct{\\uq\\h^-r))\\uo\\h^(^r), 

tG[0,T] ^ ' ^ ' ^ ' 

for some nonlinear (increasing) function Ct depending on T. 

In this paper, we begin by estimating the L^-stability for error propagation. We next 
prove that the local error of the Lie formula is an approximation of order two in time. Fi- 
nally we prove that this evolution operator represents a good approximation, of order one 



in time, of the evolution operator in the sense of Theorem 1.2 



Furthermore, we apply Lie and Strang approximations in order to make some numerical 
simulations using the split-step Fourier experiments. 



This paper is organised as follows. In the next section we give some properties related 
to the kernels G and K, and we prove two fractional Gronwall Lemmas. In Section [3] we 
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establish some estimates on X^, y*, and S*. In Section|4j we prove a local error 
estimate. Theorem 1.2 is proved in Section|5] We finally perform some numerical experi- 



ments which show that the Lie and Strang methods have a convergence rate in O (At) and 
O (At^)), respectively. 

Notations. 

- We denote by F the Fourier transform of / which is defined by: for all ^ G R, 
(1.11) J-/(0 - /(O / e-2--«/(x)dx. 

We denote by its inverse. 

- We denote by Ct (ci , C2 , • • • ) a generic constant, strictly positive, which depends on 
parameters ci, C2, • • • , and T. C is assumed to be a monotone increasing function of its 
arguments. 

2. Preliminaries 

2.1. Properties of the kernels. We begin by recalling the properties of kernels involved 
in the present analysis. 

Proposition 2.1 (Main properties of K, 11]). The kernel K satisfies: 

(1) Vi > 0, K{t, ■) G (R) and K e C°° (]0, oo[xR). 

(2) Vs, t > 0, K{s, •) * K{t, •) = K{s + t, ■). 

(3) yr > 0,3Ct > such that for all t e]0,T], \\d^K {t, ■) \\l2^r) Cxt-^^'^- 

(4) yr >0,3Ct>0 such that for all te]0,T], (i, •) IUmR-) < C*?^ 

(5) For any uq £ L^(R) and t > 0, 

\\K{t,-)*uo\\L2CR) ^e"«*||Mo|U2(R), 
where = — minRe('!/'i) > 0. 

Proposition 2.2 (Main properties of G, I9j). The kernel G satisfies: 

(!) GeC°° (]0,oo[xR). 

(2) Vs, t > 0, G(s, •) * Git, •) = G(s + t, •). 

(3) Vi>0, ||G(t,-)|Ui(R) -1. 

(4) 3Co > sMc/z that for all t > 0, 119:^6 (t, •) ||i2(R) < Got"^/''. 

(5) 3Gi > such that for all t> 0, ||a:rG(t, •) ||li(r) ^ Gii"i/2_ 

Remark 2.3. The kernel D of I — T]d^ has similar properties to the kernel K. Moreover, 
for all t > 0, we have 

(2.1) D{t,-)*G{t,-) = K{t,-). 

2.2. Fractional Gronwall lemmas. 

Lemma 2.4 (Fractional Gronwall Lemma). Let(f) : [0,T] H.^ be a bounded measurable 
function, and suppose that there are positive constants A, L and 6 Cz]0,l[ such that for all 
t£[0,T], 

(2.2) (^it)^A + L^^it), 
where is the Riemann-Liouville operator defined by 



SPLITTING FOR THE FOWLER EQUATION 



Then there exists Ct{0) such that 

< e^^('')*^, Vie[0,T]. 

Proof. The proof of this Lemma is well-known and is based on an iteration argument; see 
for instance |8, Lemma 4.3] or [22, Corollary 2]. We sketch the argument for the sake of 
completeness. Iterating inequality ( |2.2| l once, we have 

L 



^^^W)f^^' {^^W) I ^)''" V(0 dr ] ds 

= + J^tA + f\t - sf-' I {s - r)'^-V(r) drds. 



em J r(0)2 J, 

From Fubini's Theorem, we get 



{t-sf~^ / {s~rf-^(P{r)drds= / 0(r) / {t ~ sf-\s - rf-^ds dr 

Jo Jq Jr 

= 0(r)(t - r)2«-i (^J\l - rf-^r'-'dT^ dr 

= P{e,e) f <j,{r){t-rf'-Hr, 
Jo 

where (3 is the beta function. Therefore, we have 

(2.3) </)(i) ^ Ct{0)A + -^(3{0, 0) f ^{s){t - sf'-^ ds. 



n0?' 

Iterating the estimate (|2.3|l n times, with nO ^ 1, we get the following estimate: 



0(t) < Ct{.0)A + Lt{0) f ^{s){t - s)"ds, 
Jo 

with a ^ 0, and where Lt (0) isa positive constant which depends on T and 0. The lemma 
then follows from the classical Gronwall Lemma. □ 



Lemma 2.5 (Modified fractional Gronwall Lemma). Let : [0, T] — > R+ be a bounded 
measurable function and P be a polynomial with positive coefficients and no constant term. 
We assume there exists two positive constants C and 6 g]0, 1[ such that for all t e [0, T], 

(2.4) OsC0(t)<0(O) + P(i) + C7^0(O. 

Then there exists Ct{0) such that for all t G [0, T], 



m^cT{0)m+CT{0)p{t). 



6 



A. BOUHARGUANE AND R. CARLES 



Proof. Ai-guing as in the proof of Lemma 24 we iterate the previous inequality. After one 
iteration, we get 

C 







m 



\Jo 

where we have used the assumptions on P, and Fubini's Theorem again for the last term. 
Iterating sufficiently many times, we infer Uke in the proof of Lemma [24] 

(2.5) m ^ co(^(0) + coP{t) +C f (f>is){t - s^ds, 

Jo 

with a > 0. Set 



V'(i) = [co^iO) + coP{t) + CJ 0(s)(t - s)"ds ) e-^i* 

Then 

ip'{t) = (coP'(t) + Ca f (l){s){t - s)°'-^ds 



Ci (co</)(0) + c„P{t) +cj^ </)(s)(i - sr 



"ds 



Using ( |2.5| ) to control the second term, and choosing Ci sufficiently large, we infer: 

^b'it) ^ coP'{t)e-^^'. 
Since P(0) = 0, V(0) = co0(O), for all t e [0, T], 

(p{t) ^ VWe'^'* < co(/)(0)e^^* + Co / P' {s)e^'^'-'Us 

Jq 

^ co0(O)e'='i^ + coe^i^ [ P'{s) ds < coe^''^ {(j){0) + P{t)) . 
Jq 

This completes the proof. □ 



3. Estimates on the various flows 

3.1. Estimates on linear flows. In this paragraph, we collect several estimates concerning 
the convolutions with D, K and G, which will be useful in the estimates of the local error 
of the scheme. 

Proposition 3.1. Let s eILand(pe H^CR). Thenl[ip\ e H''-*^^{Il) and we have 
(3.1) ||lM|U.-4/3(R) 4^2p Q\ ii^ii^^^j^). 
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Proof. For all s £ R and all ip G H^{R), we have, using \\3) 



1/2 



3/ yjR 
2 



1 . .,.V3 
2-*sgn(0 — 



3; V-^R Vl + I^l' 
2 



\ 1/2 



(i + ieiTi-^M(ordc 



(i + ierri-^M(e)rrfe 



R 



1/2 



47r^r ( - ) II^IU.(R), 



hence the result. 



Lemma 3.2. (1) Lef n e N. Then, for all v e H"{'R) and all t > 0, 
where = — minRe(0x) > 0. 

(2) Let n S N. There exists C such that for all v G _ff^(R) and all t > 0, 
(3.2) \\X'v - i;||ff.(R) Cte^«*||«||^„+2(R). 

Proof. Using Plancherel formula, we have 

ll^*^^lli^(R) = P(i,-)*«lli2(R) 

= \\T{D{t,-))Fv\\l.^^)= I \F{D{t,-)){0\''\:Fv{0?d^ 
= [ e-^**-(«)|.F.(0|^de<e2^°lHli^(R)- 

JR 

Moreover, since 

d:x'v^D{t,-)*d^v 
then, from again Plancherel formula, we have 

hence the first point of the lemma. 
Let n e N, D e ir"+2(R). We have 

\\X'v-v\\h^^r)= f X'vds 

But from the definition of X*, X'^ is given by 

X' rjdlX'v - I[X'v] = rjX'dlv - IlX'v] 
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since X^d^v = D{s, •) * d^v — {D{s, •) * v). Thus, using Proposition 
point of this lemma, we get 



3.1 



and the first 



\X'v-v\\H„^R)^r) f \\X'dlv\\H^.^n)ds+ f \\I[X'v]\\H,.^R)ds 
Jo Jo 

<r,te'^«*||«||H"+2(R)+ / \\I[X'v]\\H^^R)ds 
Jo 

< 77te'^"*||«||H"+^(R) +4^2p ^ \\X%\\H.+2^R)ds 

< f77 + 47r2rf^))te^°*||t;||H"+2(R), 



.3^ , 

hence the result. □ 
Recalling that K corresponds to D in the case ij — 1, we readily infer: 

Corollary 3.3. For all w G iJ^(R) and all t > 0, 

(3.3) \\K{t,-)*w-w\\L2^R)^Cte''°'\\w\\H^R), 

where C is a positive constant independent oft and w. 

We conclude this paragraph with an analogous result on the heat kernel G: 

Lemma 3.4. For all w e (R) and all t > 0, 

\\G{t, ■)*w~ w||l2(r) < et\\w\\H2(^R). 

Proof. Proceeding as above, we have: 

G{t,-)*w-w = f ^{G{s,-)*w) ds = e f dl{G{s, ■) * w) ds 
Jo Jo 

= £ G{s, •) * d^w ds. 
Jo 



Taking the norm and using Proposition 2.2 Young's inequality yields 

\\G{t,-) *w- w\\l2{r) s$ e / \\G{s,-)\\LiiR)\\dlw\\L2(R) ds < et\\w\\H^R), 

Jo 

hence the result. □ 
3.2. Estimates on F*. We now turn to the viscous Burgers' equation ( |1.7| i: 
(3.4) dfW — sd'^w + wdxW = 0; W|t=o = f^o- 

Remark 3.5 (Hopf-Cole transform). It is well-known that the change of unknown function 

w = -2e\dx4> = -2edx (ln0) , 
9 

turns the viscous Burgers' equation into the heat equation: 
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We infer the explicit formula: 

However, this formula does not seem very helpful in order to establish Proposition [3]6] 

Proposition 3.6. Letwo e i/^(R). Then ( |3.4] i has a unique solution w e C(R+; if^(R)). 
In addition, there exists C — C{e, ||wo||if i(r)) such that for all t ^ 0, 

lk(i)IU2(R) ||wo||l2(r), ||5:rW(i)||L2(R) 1 1 1 1 (R) c'^'* ^ 

If in addition wq € H"'(R),for some n ^ 2, then w G C(R+; ^f"(R)) and for all T > 0, 
there exists M = M{e, T, |lwo||/f2(R)) such that for all t G [0, T], 

Proof. The existence and uniqueness part being standard, we focus on the estimates. The 
estimate yields (formally, multiply p.4| i by w and integrate) 

(3.5) ~\\Ht)\\l,+e\\dMml^=^, 

and the estimate (differentiate p.4| i with respect to x, multiply by d^w and integrate), 

(3.6) drr. 
The estimate p.5| l shows that the map 1 1 (01 1 is non-increasing: 

An integration by parts and Cauchy-Schwarz inequality then yield 

(3.7) \\dMt)\\l- < \\w{t)\\LA\dlw{t)\\L- < \\w4LA\dMi)\\L^- 

In order to take advantage of the smoothing effect provided by the viscous part, integrate 
in time and write 

' \\dlw{s)\\l.ds < \\\w',\\h + \f^ \\dM^)\\l^ds. 
Gagliardo-Nirenberg inequality yields 



so using p.7[), we infer: 



\d.M\L^^C\\dM\%^\\dlwtj^', 



In view of Holder inequality in the last integral in time, 

/■* 1 / /■* \ ''^^ 

1^ \\dMs)\\hds ^ -WuQl^ + C\\wo\\%' ^ Wdlwmhdsj t^/^ 
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where we have used Young inequality ab < a^^^ + 6®. We infer 



(3.8) 



\dM-'^)\\hd.s^\\w',\\i2+Ci\\woh^)t 



Gagliardo-Nirenberg inequahty ^ \/2||/||^2^||/'||^2'^ now yields 

\\dMs)\\L^ds < V2 / \\dMs)\\L2\\dlwis)\\%^ds 



^Cie,\\wo\\ 



L2 

t 



\dlw{s)\\%^ds 



/ft \ ''Z " 



where we have used p.7| l, Holder inequality and ( |3.8| l, successively. 

Integrate the i/^ estimate ( |3.6| l with respect to time, and now discard the viscous part 
whose contribution is non-negative: 



(3.9) 



\dxw{s)\\\-j.ds 



< Ikolli^ + 2 / \\dM^)\\LA\dxw{s)\\i,ds. 



The first part of the proposition then follows from the Gronwall lemma. 

To complete the proof of the proposition, we use the general iJ" estimate, for n G N: 
set A = (1 - 92)1/2. Applying A" to ([14]) yields 

--||A"u;(<)||i2 = 



A"w(t, a;)A" [wd^w) {t, x)dx+e / A"u;(t, x)A"5>(i, x)da;. 



R 



R 



Integrating by parts, the last term is non-positive, since 



A"w(t, a;)A"(9>(t, a;)da; = - / \K''dxw(t,x)\" dx. 



R 



R 



Write 



A"w(i, a;)A" (w^^w) (t, a;)dx = / A"w(i, x) (w^^^A'^) (t, a;)da; 



R 



R 



R 



+ / K''w{t,x){wda,K''w - K'^ {wdxw)){t,x)dx. 
J 

Integrating by parts the first term yields 



A"u;(t, x) {wd^N^w) {t, x)dx 



R 



{A"'w{t,x)) w{t,x)dx 



R 



{A''w{t,x)) dxw{t,x)dx 

R 

dxW{t)\\Loo 



In view of Kato-Ponce estimate lfT6ll 

(3.10) ||A"(/5) - /A"g|U2 < C||9./|U» + WfMah^ 
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we have (with f = w and g = dxw) 

Leaving out the viscous term, Gronwall lemma yields the a priori estimate 

(3.11) \\w{t)\\H^ i^\\wo\\w.c^v{c \\dMs)\\L^d^, 

where C depends only on n E N. In particular, for n = 2, Gronwall lemma implies 

\\d'^w{t)\\L2 < ||wo||H2e'^(*''''''+*), where C = C(e, ||wo||ffi)- We bootstrap, thanks to 
Gagliardo-Nirenberg inequality again: 

\\dMt)\\L- ^ V2\\dMt)\\L2\\dlwmL2 ^ 721^0 II ff^e^^*'''^*^. 
Therefore, f or t e [0,T], 

rt 

\\dxw{s)\\L^ds < V2\\wo\\h^ X t X e'^(^'^'+'^). 
The last estimates of the proposition then follow from p.ll| i. □ 
Lemma 3.7. Let T > 0. For all w £ i?^(R), there exists C — C{T, 111011^2 (r^)) such that 

(3.12) \\Y'w\\H^iTl) ^ e^*||ii;||ffi(K), V< G [0,r]. 
Proof. Differentiating the Duhamel formula ( |1.8[ ) in space, we have 

d^Y^w = G{t, •) *dxw- [ dxG{t - s, •) * {Y''w)dx{Y''w) ds. 



Using Young inequality and inequality ( |1.9[ ), we infer: 

\\dxY^w\\L2(Tl) l|c'a;W||L2(R) 

+ / \\d.Mt - s,-)\\LHTi)\\{Y'w)dx{Y'w)\\Li^^)ds. 
Jo 

In view of Proposition 2.2 this implies: 

\\dxY*w\\L2(R) ||<9a;W||L2(R,) 

+ Co f {t - s)-^/^\\{Y'w)dx{Y-'w)\\LiiR)ds. 
Jo 

Writing 

\\{Y''w)dx{Y''w)\\L. < \\Y''w\\L2\\dx{Y'w)\\L^ ^ \\w\\l4MY'w)\\l2, 



and invoking the fractional Gronwall Lemma 2.4 with 6* = 1/4, the lemma follows. □ 
Corollary 3.8. Let n^landw £ iJ"(R). Let T,a>0. If 

I|w||l2(r) ^ a, 
then there exists c depending only on T and a such that 

\\Y'w\\h^.^r) e^*||u;||H,.(R), Vt G [0,T]. 
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Proof. Denote w(t) = Y^w. From Lemma 



3.7 



^(<)l|ffi(R) < ||wo||j/i 



where C depends only on a and T. From the proof of the first part of Proposition 3.6 (see 
the estimate after p.8| l), we infer 

t 

\\dMs)\\L^ds^C{a,T), Vt e [0,T]. 
The corollary then stems from ( |3.11| l. □ 

3.3. Estimates on the splitting operator Zj^. Combining the estimates on X* and F* 
established in the previous two sections, we infer: 

Corollary 3.9. (1) For all u e L^(R) and all t > 0, 

\\ZIu\\l2i^r) < e'3''*||u||i2(R), 

where = — mmKe{(f>x) > 0. 
(2) Let r > 0, n € N* and u € i/"(R). There exists C = C{T, ||u||l2(r)) such that 
for all t £ [0, T], 

II^l"IIh"(R) «^ e'^*ll"llff"(R)- 



Proof. The first point is a direct consequence of the relation (|1.9|l and Lemma 3.2 



The second point is readily established with Lemma 3.2 and Corollary 3.8 □ 
3.4. Estimates on the exact flow S*. 

Lemma 3.10 (i^-a priori estimate). Let uq S ^^(R.) and T > 0. Then, the unique mild 
solution u e C([0,r];L2(R)) n C(]0, T]; i72(R)) of (pTji satisfies, for all t £ [0,T] 

\\u{t, •)IU2(R) sC e"»*||wo||L2(R), 
where = — minRe('0x) > 0. 

Proof. Multiplying \\.\\ by u and integrating with respect to the space variable, we get: 



UfU dx + / {I[u\ ~ Uxx) u dx = 

R Jr 

because the nonlinear term is zero. Using ( |1.3| l and the fact that u and J-^{I[u] — dl^u)u dx 
are real, we get 

/ {I[u]^dl^u)udx^ I F-^{i;iFu)udx [ ipxlJ'uf d^ = [ Re{^|Jx)\J'uf d£,. 
Jr Jr Jr Jr 

We infer 

where uq — — min Rc(V'i) > 0. The result then follows from the Gronwall lemma. □ 

Lemma 3.11. Let n £ N*, uq e i7"(R) and T > 0. There exists CT(||uo||ff"-i(R)) 
such that the unique mild solution u e C([0, T]; iJ"(R)) satisfies 

(3.13) \\u{t,-)\\H'-{R) < C'T(||wo||if"-i(R))||wo||H"(R)- 
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Proof. The proof is similar to the one given in Lemma 3.7 Differentiating the Duhamel 
formula \ \A\ in space, we have 

d^u = K{t, •) * dxUo - / dxK{t - s, •) * (u(s)9^u(s)) ds. 



Using Young inequality and Proposition |2T| we infer, for any integer n ^ 1: 

\\dxUit)\\H^-i{-R.) < e"°*||a:rMol!ff"-i(R) 



In view of Proposition 2. 1 this implies: 

\\dMt)\\H'^-HR) ^ e"°*||c)^uollH"-i(R) 



Jo 

For n = 1, we use Lemma |3.10| to have 

\\u{s)dxu{s)\\Li < \\u{s)\\L2\\d.j;u{s)\\L2 < e'^"'^\\uQ\\L2\\dxu{s)\\L2. 



The fractional Gronwall Lemma 2.4 with = 1/4 then yields 



where C depends only on T and ||uo||l2. From ( |L9[ ), this implies the lemma in the case 
n = 1. For n ^ 1, Leibniz rule and Cauchy-Schwarz inequality yield 

||m(s)9j;U(s)||vi/"-1-1(R) < C{n)\\u{s)\\H'^-i(R)\\dxU{s)\\H^-i(Ry 

The lemma then easily follows by induction on n. □ 
We will also need the fact that the flow map S* is uniformly Lipschitzean on balls of 

i?2(R). 

Proposition 3.12. Let T,R> 0. There exists K — K{R, T) < oo such that if 

\\uo\\m{R) ^ R, \\vo\\hHR) < R, 

then 

\\S*'uo - S'*vo||l2(r) ^ K\\ua - wo||l2(r), Vt G [0,r]. 
Proof. Set u{t) = S^uq, v{t) = S^vq and w = u — v. It solves 
(3.14) dfW +I[w] — d^^w = vdxV — udxU — —udxW — wdxV. 

The energy estimate yields: 



\ d f 

2^II'^IIl2+ J w{udxW + wdxv) ao\\w{t)\\'l2, 



We 



where the term J w {l[w] — d^^w) has been estimated as in the proof of Lemma 3.10 
infer 

^^ii^i'Wiii^ < («o + l\\dxu{t)\\L^ + \\dxvit)\\L^^ Wwmh 

^ c(i + \\umH2 + ikwiko ii^wiii^ cxR,T)\\wmh, 

where we have used Sobolev embedding and Lemma |3.11| Gronwall lemma yields the 
result. □ 
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4. Local error estimate 



As pointed out in the introduction, the formalism of Lie derivatives does not seem to 
be easy to use in the framework of this paper, since the action of the nonlocal operator I 
on nonlinearities is rather involved. As a consequence, we prove an error estimate by a 
rather pedestrian way. 

Proposition 4.1 (L^ local error estimate). Let uq e if^(R). There exists C (\\uo\\l^(R)) 
such that for all t € [0, 1], 

\\Zluo - S'*uoI1l2(r) C (||moIIl2(r)) \\uo\\h^r) 



Proof. From the definition of Z]^ and Remark 



2.3 



!■)■ 

we have 

t 

2 



Git~s)*d^{Y'uoYds 



1 



= D{t) * G{t) * uo - ^ / D{t) * Git - s) * d^iVuQ)'^ ds 



1 



(4.1) = K(t) *uo-^ I D{t) * Git - s) * d^iY'uof ds. 



Thus, from Duhamel formula for the Fowler equation ( |1.4| i and the Lie formula ( |4.1| i, we 
have: 

Zluo - S'uo = 1 f d^Kit - s) * iS'uof ds- \ j Dit) * O^Git - s) * iV'uof ds 
^ Jo ^ Jo 



1 



(4.2) =^j^ d^Kit - s) * {iS'u^f - iZiuof) ds + Rit), 

where the remainder is written as 



t 

2 



Rit) = I I Riis)ds, with = d,Kit-s)^iZluoy-Dit)^d^Git~s, XF^uq)' 

^ Jo 

Then, from Proposition 2.1 Corollary 3.9 and Lemma [3.10| we have, for t e [0, 1]: 



1 



\\Zluo~S'uo\\mii) < 27 WdxKit- s,-)\\mn)\\iS'uor - iZluor\\LHTL)ds 
+ \\Rit)\\mR) 

f it- s)-3/4||5-yQ _ ZIuo\\l2^r)\\S'uo + ZIuoWlhr) ds + \\Rit)\\mu.) 
Jo 

<: c(e"°* + e^«*)||uo||L^(R) f it- sr^'^^S'uo - Zluoh-^R) ds + \\Rit)\\mR), 



where C is a positive constant. To estimate the remainder, we decompose it as follows 

i?i(s) = Ti+T2+T3+Ti, 

where 

Ti = Kit - s, •) * d^iZluof - d^iZiuof, 
T2 = d^iY'uof - Git -s,-)* d^iY'uof, 
T3 = Git ~s,-)* dAY'uof - Dit, •) * Git - s, •) * d^iY-'uof, 
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Let us first study the term Ti. From Corollaries 3.3 and 3.9 we have (recall that t E [0, 1]) 

I|Ti||l^(r) = \\K{t - s, •) * d^iZluof - a,(Z£wo)'||L2(R) 

<Ce"°(*-^)(t-s)||Zi^.o||?f3(R) 
<C(ho||L^(R))e"°(*-^)(i-s)||zio|ll^3(R). 
In the same way, from Lemma [X?] and Corollary |3. 8 1 we control the term T2 as 

||r2|U2(R) = Wd^iY'uof - G{t - s, •) * d^Y^uofh^t^R) 

L2(R)) {t - s)||uo||^3(R)- 

From Lemma [J!2l and Corollary |3.8| 

mWmR) = \\G{t - •) * d^iY'uo f - D{t, •) * G{t -s,-)* d,{Y'uonmR) 
Ce^o' t \\G{t - s, •) * dAY'uofWHHR) 
^Ce^'>U\\d.,{Y'uof\\HHii) 

ff3(R)- 

For the term T4, write 

I|74||l2(r) = 119^ {Zluof - (F'^uo)^ IU2(R) 

= 2||(Z£uo)9x (^£t^o) - (i"^Ho)9. (y^«o) ||l2(r). 
By linearity of the evolution operator X*', we have 

a, (z£mo) = ^'9, (r^uo) , 

hence 

I|74||l2(r) - 2||(Z£uo)X^a, (r^uo) - (l^^wo)9. (r^uo) IIl-(r) 
^ 2 11X^9, (y^uo) [X'Y'uQ - >''^o)|Il2(r) 

+ 2 ||(r-^«o) (^^a, (r^uo) - 9, (r^uo))lli2(R) . 

Now from Sobolev embedding. Lemma [l!2] and Corollary |3.8| we get: 

\\Ti\\LHTl) < 2||X^9, {Y'uo) ||L-(R)||X^r^Uo - r^^^o||L2(R) 

+ 2||r^uo||L-(R)ll^'5, (y'«o) - d., {Y'uo) |U2(R) 
^ C\\X'd^ (Y'uq) |Ui(r) e'^"^- s \\Y'u^\\m(Yi) 
+ C||r^«o||Hi(R) e^«^ 5 ||9, (r^uo) IIh2(r) 

<C(||zio||L2(R))e2''"^s||uo|Pff3(R). 
Finally, since i?i (s) = Ti + T2 + Tg + T4 then for < s < 1, 

Pi(s)IIl2(r) < C(||uo||l2(r)) i ||uo|Ih3(r), 
and by integration for s e [0, t], 

ll-RWIU2(R) < (|1uo||l2(r)) ||uo||ff3(R). 
We conclude by applying the modified fractional Gronwall Lemma 2.5 □ 
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5. Proof OF Theorem [TT2 



The proof follows the same lines as in ifTSl Section 5]. Denote by = (-^l *) '^o the 
numerical solution, and 

which corresponds to the exact evolution of the numerical value Uk at time tk — kAt up 
to time tn = nAt. From Lemma [3.10| there exists p such that 

II«WIIh^(r) s^p, vte[o,r]. 

We prove by induction that there exists 7, Ato, c > such that if < At ^ Ato, for all 
n G N with nAt ^ T, 

\\un\\m(R) ^ e™'^*||uoll_f/3(R) < Co, 

where Co = e'^-'"||uo||_ff3. The above properties are satisfied for n ~ 0. Let n ^ 1, and 
suppose that the induction assumption is true for ^ k ^ n — 1. Since u„ = and 
u{tn) = w", we estimate 



n~l 



\\Un - u{tn)\\L2 ^ ^ II 



k+1 



k=0 
n-1 



k=0 



L2 



For fc ^ n — 2, Z^^Uk — Ufc+i and Proposition 3.12 yields, along with the induction 
assumption, 

^ ifllufc - u(tfc)||L2 + ||u(tfe+i)||i2 K-iAt + p, 

which is bounded by 2p if < At ^ Ato ^ 1- Up to replacing K with max(i4r, 1), we 
obtain, for /c ^ n — 1 and nAt ^ T, 



Using Proposition 4. 1 we infer 



^(„-fc-i)At {^zt'uk) - , ^ Cif(At)2e2='=^*||uo||ff3, 



for some uniform constant C. Therefore, 

\\un - u(t„)||L2 sC nCA'(At)2e2'='^||uo||lj3 s$ CTKe^^^At, 

which yields the first two estimates of the induction, provided one takes 7 CTKe^'^'^, 
which is uniform in n and At. Finally, the last estimate of the induction stems from 
Corollary 
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6. Numerical Experiments 

The aim of this section is to numerically verify the Lie method convergence rate in 
O (At) for the Fowler equation 

To solve the linear sub-equation ( |1.6| l, discrete Fourier transform is used and for the non- 
linear sub-equation ( |1.7| i, different numerical approximations can be used. Here, we use 
the finite difference method. 

Since the discrete Fourier transform plays a key role in these schemes, we briefly review 
its definition, which can be found in most books. In some situation, when the mesh nodes 
number N is chosen to be = 2^ for some integer p, a fast Fourier transform (FFT) 
algorithm is used to further decrease the computation time. In this work we will use a 
subroutine implemented in Matlab. In this program, the interval [0, 1] is discretized by N 
equidistant points, with spacing Aa; = 1/iV. The spatial grid points are then given by 
Xj — j/N, j — 0, N. If Uj{t) denotes the approximate solution to u{t, xj), the discrete 
Fourier transform of the sequence {uj }^Sq is defined by 



for fc = 0, • • • , — 1, and the inverse discrete Fourier transform is given by 



1 



N 

k=0 



for j = 0, • • • , TV — 1. Here denotes the discrete Fourier transform and T its inverse. 



In what follows, the linear equation ( |1.6[ ) is solved using the discrete Fourier transform 
and time marching is performed exactly according to 

(6.1) = F-"^ {e-^^^''^^*'Ft{u' 



To approximate the viscous Burgers' equation we use the following explicit centered 
scheme: 



(6.2) 



^ ' 2Ax 



2 / V 2 , 



+ eAt 



which is stable under the CFL-Peclet condition 

'Ax Ax^' 

where v is an average value of u in the neighbourhood of (t", Xj). 



(6.3) At = min 



Remark 6.1. In the case where the linear sub-equation ( |1.6| l is solved using a finite differ- 
ence scheme instead of a FFT computation, an additional stability condition is required, 
see Is). Moreover, the computation time becomes very long because of the discretiza- 
tion of the nonlocal term which is approximated using a quadrature rule. Indeed, in |3|, 
the Fowler equation has been discretized using finite difference method and the numerical 
analysis showed that this operation is computationally expensive. This observation has 
also motivated the use of splitting methods, in particular the implementation of split-step 
Fourier methods. 
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In order to avoid numerical reflections due to boundaries conditions and to justify the 
use of the FFT method, we consider initial data with compact support displayed in Figure[T] 
to perform numerical simulations. 




-1 -Z -, . , Z 1 -1 -z -1 1 z 1 -1 -z -1 1 z 

(a) (b) (=> 

Figure 1 . Initial data used for numerical experiments. 

Since we do not know the exact solution of the Fowler equation, a classical numerical 
way to determine the convergence numerical order of schemes is to plot the logarithm of 
the error ||wi(r) — U2(T)||i2 in function of the logarithm of the step time At, where ui 
and U2 are computed for time steps At/ 2 and At /A, respectively, up to the final time T. 
Hence, the numerical order corresponds to the slope of the curve, see Figures |2] |3] For 
reference, a small line of slope one (resp. two) is added in Figure |2] (resp. [3]l. We see 
that the slopes for the three initial data match well and so we can conclude that numerical 
simulations are consistent with the theoretical results estabUshed above. 




log(At) 

Figure 2. Lie method 
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We also study numerical convergence of Strang splittings using initial data displayed in 
Figure[T| Results are plotted in Figure|3] We can see that the Strang formulation is of order 
two in time for smooth initial data. 




log(At) 

Figures. Strang method 

From numerical simulations, we emphasize the fact that both formulas defining a Lie 
operator, as well as both formulas defining a Strang operator, lead to the same results. 

Acknowledgements. The authors are grateful to Pascal Azerad and Bijan Mohammadi 
for helpful comments. 
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